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Optical Feshbach resonances [Phys. Rev. Lett. 94, 193001 (2005)] and pump-dump photoas- 
sociation with short laser pulses [Phys. Rev. A 73, 033408 (2006)] have been proposed as means 
to coherently form stable ultracold alkali dimer molecules. In an optical Feshbach resonance, the 
intensity and possibly frequency of a cw laser are ramped up linearly followed by a sudden switch-off 
of the laser. This is applicable to tightly trapped atom pairs. In short-pulse photoassociation, the 
pump pulse forms a wave-packet in an electronically excited state. The ensuing dynamics carry the 
wave-packet to shorter internuclear distances where, after half a vibrational period, it can be deex- 
cited to the electronic ground state by the dump pulse. Short-pulse photoassociation is suited for 
both shallow and tight traps. The applicability of these two means to produce ultracold molecules 
is investigated here for ***Sr. Dipole-allowed transitions proceeding via the B^T,i excited state as 
well as transitions near the intercombination line are studied. 

PACS numbers: 33.80.-b,32.80.Qk,34.50.Rk,33.90.-l-h 
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I. INTRODUCTION 

The cooling and trapping of alkaline-earth metals and 
systems with similar electronic structure such as ytter- 
bium have been the subject of intense research over the 
last decade. The interest in ultracold group II atoms was 
triggered by the quest for new optical frequency stan- 
dards [H. The extremely narrow linewidth of the inter- 
combination transition, together with the magic wave- 
length of an optical lattice , is at the heart of the clock 
proposals. This narrow width has several interesting im- 
plications, such as a very low Doppler temperature for 
laser cooling @, or optical Feshbach resonances involv- 
ing small losses and thus holding the promise of easy 
optical control [3, 's', '3] . Additionally ultracold alkaline- 
earth metals offer many more exciting perspectives: They 
are candidates for high precision measurements 0, Q , for 
studies of ultracold mixtures @, for the realization of 
quantum information processing [Tol . [Tl| and quantum 
phases with unusual symmetries jl2| . 

The present study focusses on strontium and addresses 
the question of how to produce dimer molecules in their 
electronic ground state from ultracold ®^Sr atoms. Em- 
phasis is put on a coherent formation mechanism. Two 
different time-dependent schemes are studied, a ramp 
over an optical Feshbach resonance and a sequence 
of short pump-dump laser pulses cf. Fig.[Tl that are 
both based on photoassociation resonances. 

Photoassociation, the excitation of two colliding atoms 
into a weakly bound molecular level of an electronically 
excited state by laser light, has emerged as a standard 
tool to study the scattering properties of ultracold atoms 
[isj . For alkaline-earth species, both dipole allowed tran- 
sitions and transitions near the intercombination line can 



be employed. Photoassociation based on dipole-allowed 
transitions near the ^Sq — ^Pi atomic resonance was ob- 
served for calcium , for strontium [13, > and for yt- 
terbium |.19i] . Near the narrow-line atomic intercombina- 
tion transition, photoassociation spectroscopy has been 
performed for strontium [2l| and for ytterbium [l^j • 

In addition to gathering spectroscopic information, 
photoassociation serves as a method to produce stable 
ultracold molecules [1^. Depending on the properties of 
the electronically excited state, spontaneous decay may 
yield molecules in their electronic ground state [2^ . How- 
ever, spontaneous decay toward molecules instead of a 
pair of atoms occurs only if the vibrational wavefunc- 
tions show a large probability amplitude at short in- 
ternuclear distances R. On the other hand, large free- 
bound Franck-Condon factors and hence very extended 
vibrational wavefunctions are required for photoassocia- 
tion. Ideally the vibrational wavefunctions should there- 
fore consist of both a short-range part and a long-range 
part. In heavy alkali dimers such wavefunctions exist 
thanks to strong spin-orbit coupling ^23]. In particular. 
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FIG. 1: (Color online) Coherent formation of ultracold stron- 
tium molecules in their electronic ground state via (a) an 
optical Feshbach resonance where an atom pair in the lowest 
state of a tight trap is adiabatically compressed to a molecule 
and (b) by short pulse photoassociation where molecules are 
formed in a pump-dump sequence. 
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the softly repulsive part of a purely lon g-ra nge well (23| 
and resonant spin-orbit coupling [3) Izsl l26j can provide 
the necessary 'i?-transfer' mechanism. If the vibrational 
wavefunctions of the electronically excited state do not 
possess this dual character, molecules in the last one or 
two levels of the electronic ground state can be formed 
at most. Overwhelmingly, spontaneous decay will simply 
redissociate the molecules to a pair of atoms. 

Photoassociation as a means to produce molecules in 
their electronic ground state may be more efficient if 
time-dependent laser fields are employed. For example 
a sequence of two laser pulses, a pump pulse and a time- 
delayed dump pulse can be utihzed (see also Fig.[T]D). 
The pump pulse excites a vibrational wave-packet in the 
electronically excited state. The wave-packet moves to 
shorter internuclear distances where it can be transferred 
to the electronic ground state by the dump pulse. The 
character of the vibrational wavefunctions decides the 
fate of molecule formation: The wave-packet dynamics 
lead to appreciable overlap with molecular levels in the 
electronic ground state only if the eigenfunctions from 
which the wave-packet is built provide for it [27]. This 
reflects the weak field approach implicit in the pump- 
dump scheme where the maximum transition probabil- 
ities are completely determined by the Franck-Condon 
factors and the spectral amplitudes of the field. 

An alternative time-dependent scheme for creating 
molecules in their electronic ground state consists in 
ramping the intensity and possibly also frequency of a 
continuous-wave (cw) laser over a photoassociation res- 
onance [i^ (see also Fig. [T^). This has been termed 
optical Feshbach resonance (FR), in analogy to the for- 
mation of molecules with magnetic fields 28] . It requires 
an initial condition where two atoms are not too far from 
each other such as in a Mott insulator state of an optical 
lattice with two atoms per lattice site. The intensity of 
the laser field is ramped up slowly while the frequency 
of the field is chosen such that the population of the 
excited state is kept minimal 13|. This ensures an adia- 
batic compression of the wavefunction while avoiding loss 
due to spontaneous emission. In the case of rubidium, a 
molecule formation efficiency of 50% was predicted for 
realistic laser parameters. 

While the pump-dump sequence with short pulses and 
the ramp over an optical FR with a cw laser both rep- 
resent coherent molecule formation schemes, they differ 
in a number of aspects. (1) The timescale for the pump- 
dump sequence is determined by the vibrational periods 
in the electronically excited state which are on the or- 
der of 100 ps. Ramping over an optical FR needs to be 
done much more slowly, typically over 0.1 — 1 /iS. This 
timescale is set by the initial state, i.e. by the frequency 
of the trapping potential. (2) Short pulses address only 
a small part of the initial wavefunction in a certain range 
of internuclear distances, the photoassociation window. 



When ramping over an optical FR, the whole wavefunc- 
tion is adiabatically deformed. 

Here, the two schemes are applied to the formation of 
strontium dimers in their electronic ground state. They 
are schematically depicted for optically allowed transi- 
tions proceeding via the Sr2 -B^Sj excited state in Fig.[T] 
Whether coherent molecule formation by optical transi- 
tions near the intercombination is possible, is also stud- 
ied. Similarities and differences between alkali and al- 
kaline earth dimers are analyzed by comparing *^Sr2 to 
^^Rb2. The paper is organized as follows. Section [II] 
outlines the two-state model that is required to describe 
molecule formation by optically allowed transitions as 
well as the three-state model that presents a simplified, 
yet sufficiently accurate view of the physics near the in- 
tercombination line. The optical FRs are studied in Sec- 
tion mil short-pulse PA is discussed in Section IIVI and 
Section IV] conludes. 



II. MODEL 

Two strontium atoms initially interacting via the elec- 
tronic ground state potential are considered. Assuming 
a harmonic trapping potential, the center of mass mo- 
tion can be integrated analytically and is omitted from 
the description. For sufficiently low temperatures only 
s-wave encounters need to be taken into account. The 
electronic ground state potential supports both bound 
levels and scattering states. A laser can excite an un- 
bound atom pair via an electric dipole-allowed transition 
into bound molecular levels below the ^S'o -Pi asymp- 
tote, or via transitions near the intercombination line into 
bound molecular levels below the ^Sq +^ Pi asymptote. 

The Hamiltonian modelling electric dipole-allowed 
transitions between the A"^S^ electronic ground state 
and the -B^S+ electronically excited state is given by 



ALit) ~ f F 



(1) 



where the rotating-wave approximation (RWA) is in- 
voked. The single-channel Hamiltonians, H£,(e) , consist of 
the kinetic energy, the molecular potential energy curves 
and the trap potential, Hg(e) = T-l- Vg(e) (R) -I- (— )Vtr(R). 
The potentials Vg(e) are schematically shown in Fig. (TJ 
they are taken from Refs. [Hj (Vg) and [H (Ve). The 

trap is assumed to be harmonic, Vtr(R) = ^muj^^R with 
TO denoting the reduced mass and Wtr the trapping fre- 
quency. 

Following Ref . [13] , a three-channel Hamiltonian, again 
in the RWA, is used to model transitions near the inter- 
combination line, 
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AL(i)-fr, 



(2) 





near the intercombination line 


dipole allowed 


C3 

lifetime 


0.0075 Eh 
21.46 Ais 


18.54 Eh al 
5.22 ns 



TABLE I: C3 coefficient and lifetime of strontium atoms for 
electronic states accessed by transitions near the intercombi- 
nation line and by dipole-allowed transitions 



with the rotational couplings given in terms of K-ot(R) = 

- 2 

^ . The potentials of Lennard- Jones type fitted to ex- 



perimental data are found in Ref. [20|. 

The transition dipole operators are assumed to be in- 
dependent of R and fixed at their atomic values, i.e. fi 
and p, are given in terms of the respective C3 coefficients, 
cf. Table HI The electric field amplitude is denoted by 
E{t). In case of the optical FR it corresponds to a lin- 
ear function of time or to a constant, while a Gaussian 
pulse envelope is assumed for the short-pulse photoasso- 
ciation. The detuning, Ai(i) = h{u!a,t — a;i(t)), becomes 
time-dependent only in the case of the optical FR, and 
for simplicity, a linear time-dependence of Ai(t) is cho- 
sen. The finite lifetime of the excited states is modelled 
by an exponential decay with decay rates F, F. This 
corresponds to the assumption that every excited state 
molecule which undergoes spontaneous emission yields a 
hot atom pair which is lost from the trap. The lifetimes 
corresponding to F, F are listed in Table HI 

The Hamiltonians are each represented on a Fourier 
grid employing an adaptive grid step [1^ [13, HH . Set- 
ting E{t) = 0, the wavcfunctions and binding energies 
of the vibrational levels are obtained by diagonalization. 
Analogously, dressed states are calculated by diagonal- 
ization for fixed values of E{t) and AL(f). 

Coherent optical molecule formation is investigated by 
solving the time-dependent Schrodinger equation. 
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with a Chebychev propagator 32] where H(i) corre- 
sponds to either Eq. ^ or to Eq. The state of the 
system at time t, is represented on the mapped 

Fourier grid for each channel i, 'i>i{R;t) — {R\'^i{t)). 



III. STRONTIUM MOLECULE FORMATION 
VIA OPTICAL FESHBACH RESONANCE 

Molecule formation via an optical FR is based on 
the formal analogy between a photoassociation resonance 



and a magnetic FR. Instead of ramping the magnetic field 
magnitude, the intensity and detuning of a cw-laser are 
increased in order to cross the resonance il3)]. Then, in 
the dressed-state picture, the number of eigenstates be- 
low the ground state dissociation limit is increased by 
one. If the ramp is performed slowly enough, the low- 
est trap state is adiabatically connected to the dressed 
state just below the ground state dissociation limit. This 
dressed state is an eigenstate of the Hamiltonian for the 
final values of field amplitude and detuning. However, 
the excited state components of the state of the system 
are subject to decay due to spontaneous emission. There- 
fore the laser coupling between the ground and excited 
state should be switched off before significant sponta- 
neous emission losses occur. A sudden switch-off of the 
coupling corresponds to a projection of the dressed state 
onto the eigenstates of the field-free Hamiltonian, i.e. 
onto the vibrational levels. If the dressed state before 
the switch-off has a large overlap with the last bound 
level of the electronic ground state potential, molecules 
are formed. 

The success of this scheme is determined by the value 
of the dressed-state projections onto the last bound level 
which depends on the laser intensity and detuning, and 
by the time within which an adiabatic ramp can be per- 
formed. This time is constrained from above and be- 
low by the requirement of avoiding spontaneous emission 
loss and by the requirement of adiabaticity, respectively. 
Both limiting timescales are controlled by experimentally 
accessible parameters. The loss due to spontaneous emis- 
sion depends on how much population is actually trans- 
ferred to the electronically excited state - this is deter- 
mined by the laser intensity and detuning. On the other 
hand, adiabaticity requires the ramp to be slower than 
the slowest internuclear dynamics, that is slower than the 
vibrational period of the lowest trap state ~ this is de- 
termined by the trap frequency, i.e. the intensity of the 
trapping laser. Vibrational periods in the sub-/is regime 
are realized in tight traps such as a deep optical lattice. 



A. Dipole allowed transitions X^T,g 



The two quantities characterizing gain and loss in 
molecule formation via an optical FR, the projection 
of the dressed state onto the last bound level and the 
amount of excited state population in the dressed state, 
are shown in Fig. [5] as a function of laser intensity and 
detuning. When the laser detuning coincides with the 
binding energy of a vibrational level of the electronically 
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FIG. 2: (Color online) Excited state population (bottom) and 
projection onto the last bound level (top) as a function of laser 
intensity and detuning for three different pairs of photoasso- 
ciation resonances 



excited state, large values of the projection are observed 
even at small laser intensity. Such a resonant transi- 
tion comes at the expense of transferring a significant 
amount of population to the electronically excited state. 
Rather, detunings inbetween two resonances can be ex- 
pected to yield large projections while keeping the pop- 
ulation transfer to the excited state small. Note that 
the scale of the excited state population is cut-off at 
0.01: The dark-red regions of the excited state population 
should be avoided since spontaneous emission losses will 
be too large. Three pairs of photoassociation resonances 
are studied: (i) 3.938 cm~^ and 4.135 cm~^ (left column 
ofFig.[2]), (ii) 5.006 cm~^ and 5.247 cm"^ (middle column 
of Fig. [5]), (iii) 6.306 cm^^ and 6.597 cm^^ (right column 
of Fig. [5]). As the vibrational level spacing increases from 
left to right in Fig. [21 the areas corresponding to the de- 
sired large projections and the undesired large excited 
state populations shrink. Small regions exist at large in- 
tensity and close to the photoassociation resonance where 
the dressed state projection onto the last bound level is 
large while the excited state population of the dressed 
state is small. If the final values of field intensity and 
detuning correspond to such a region, the molecule for- 
mation scheme is expected to be successful. In order to 
determine the necessary ramps in intensity and detun- 
ing, a path needs to be identified in Fig. [2] which starts 
at zero intensity and leads to the target values of field 
intensity and detuning without passing through any of 
the dark-red regions with large excited state population. 
A possible choice, studied below in Fig. [31 consists in 
ramping the intensity from zero to 1000 W/cm^ at a de- 
tuning of — 6.36cm^^, and then ramping the detuning to 
— 6.34cm~^ (corresponding to an overall change of de- 
tuning by 600 MHz). The projection onto the last bound 
level of the dressed state for the final values of inten- 
sity and detuning amounts to 51%. Larger projections 
would require a larger final intensity, and most impor- 
tantly, ramping the detuning over about 2 GHz. 

Whether molecule formation in this way is feasible re- 
lies on the existence of a window of time inbetween the 




time [ ns ] 

FIG. 3: (Color online) Projection of the time-dependent wave- 
function, |^(t)), onto the last bound level and the first trap 
states as a function of time for a ramp of intensity followed 
by a ramp of frequency (bottom panel). The field amplitude 
(solid line) and detuning (dashed line) are shown in the upper 
panel. 

timescales set by the adiabaticity condition and by the 
requirement to avoid spontaneous emission. An estimate 
of the timescale for spontaneous emission loss is obtained 
by dividing the lifetime by the maximum amount of ex- 
cited state population. For laser intensities and detun- 
ings where significant overlap of the dressed state with 
the last bound level is achieved, the excited state popu- 
lation is typically at least 1 %. Given the lifetime of Sr2 
which is by a factor of 5 shorter than that of Rb2 , an up- 
per limit in time of 370 ns is obtained, compared to 1.9 /is 
for Rb2 [l^. This needs to be related to the vibrational 
period of the lowest trap state which is 29 ns (51 ns) for 
i/t, = 500 kHz (300 kHz). 

Fig. [3] displays the total population together with the 
projections of the time-dependent wavefunction onto the 
last bound level and the first trap states as a function 
of time for the ramp described above. The overall ramp 
time is chosen to be 60 ns. Clearly, this time is too short 
for the adiabaticity condition to be fully met. However, 
in that time, already 30% of the population is lost due to 
spontaneous emission (dotted line in Fig. [3]). The losses 
occur in particular during the frequency ramp, and non- 
adiabatic oscillations in the time-dependent projection 
onto the last bound level as well as repopulation of the 
lowest trap state (red solid line and black dashed line in 
Fig. [3]) are observed. 

Ramps in intensity only are therefore studied in Fig.[4l 
The detuning is kept fixed and corresponds to that of 
the final detuning in Fig. [3] while the intensity is again 
ramped from zero to 1000 W/cm^. Comparing the left- 
hand side of Fig. [4] to Fig. [31 the population loss is re- 
duced to 20% and no indication for deviations from adi- 
abatic following are observed in the projection onto the 
last bound level (red solid line) and the lowest trap (black 
dashed line). The probability for forming molecules 
amounts to 25% in Fig. [H compared to 23% in Fig. [31 
Higher probabilities require larger final intensities. For 
example, a molecule formation probability of 33% is 
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FIG. 4: (Color online) The same as in Fig. |3]for a ramp of 
intensity only, the detuning is kept constant at — 6.36cm~^. 
Two different trap frequencies are compared with the adia- 
baticity constraint better fulfilled for the tighter trap. 



achieved for a final intensity of 1500 W/cm^ (data not 
shown). In its left-hand and right-hand sides, Fig. 0] 
compares identical ramps for different trap frequencies. 
Adiabatic following is better enforced in the tighter trap 
corresponding to the faster dynamical timescale. This 
shows up in the smaller projection onto the last bound 
level and slower depletion of the lowest trap state for 
j^tr — 300 kHz as compared to — 500 kHz . It is also 
indicated by the smaller spontaneous emission loss corre- 
sponding to less population transfer to the excited state. 

Compared to Rb2, it is more difficult to produce 
ground state molecules via an optical FR on a dipole- 
allowed transition in Sr2. The shorter excited state life- 
time of strontium requires faster ramps and leads to 
stronger loss of population. Only for very tight traps, 
adiabatic following can be realized. 



B. Transitions near the intercombination line 

Molecule formation via an optical FR near the inter- 
combination line turns out to be impossible. In the inter- 
pretation given above, an optical FR corresponds to an 
increase of the number of levels below the ground state 
dissociation limit in the dressed state picture. Generally, 
however, an optical FR occurs when the number of levels 
below the ground state dissociation limit in the dressed 
state picture is changed by one. This number of levels 
could also be decreased. In that case, the usual signature 
of a FR, a change in the scattering length, is observed. 
The dressed state wavefunction is, however, pushed out- 
ward to larger distances instead of being compressed 
inward. When crossing the resonance, the last bound 
molecular level becomes unbound, and obviously it is 
not possible to form molecules. This phenomenon cor- 
responds to open-channel dominated vs closed-channel 
dominated magnetic FRs. 

In the optical case, which direction a FR takes is deter- 
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FIG. 5: (Color online) Light shift; Binding energy of the least 
bound dressed state levels as a function of intensity for the 
laser tuned close to the second to last excited state level H 
(last level - black circles, second to last level - red squares, 
third to last level - blue diamonds, fourth to last level - green 
triangles). The upper panel shows the behavior of the last 
level near zero laser intensity and near the resonance. 



mined by the direction of the light shift. Typically, a neg- 
ative light shift corresponding to an increase in the num- 
ber of bound levels is observed for 1 /R^ long-range po- 
tentials while long-range potentials yield a positive 
light shift. The potentials accessed by transitions near 
the intercombination line show a long-range be- 

havior but with a very small C3 coefficient (cf. Table 
This leads to a curious cross-over from from dipole-dipole 
to van der Waals interaction, i.e. from 1/R'^ to 1/R^ be- 
havior (cf. Ref. [205 where the nine least-bound excited 
state levels were measured). Only the binding energies of 
the last two levels correspond to the energy range where 
the 1/R'^ terms of the excited state potentials dominate 
the 1/R^ terms. 

The binding energy of the dressed state levels, i.e. the 
light shift, is shown as a function of intensity for the laser 
tuned close to the second and fourth to last levels of the 
excited states in Figs. [5] and [6l The corresponding bind- 
ing energies are 24.21 MHz and 353.28 MHz. Note that 
the binding energies are shifted by the laser detunings, 
Al = 24.0 MHz in Fig. O and = 353.0 MHz in Fig. El 
i.e. the dressed state representation is employed also at 
zero laser intensity. The binding energies appear there- 
fore as at zero For both photoassociation resonances, the 
number of dressed state levels below the ground state 
dissociation limit is decreased rather than increased: As 
the resonance is crossed (at about 0.6kW/cm^ in Fig. O 
and at about 160kW/cm^ in Fig. [S]), the black circles 
are pushed above the dissociation limit, and the second 
to last becomes the last level: The light shift is positive 
for both resonances. The 1/i?^ character of the excited 
state potentials shows up only at very low intensity for 
the resonance at 24.21 MHz, cf. the inset of Fig. [S] The 
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FIG. 6: (Color online) Same as Fig. [5] but for the laser tuned 
close to the fourth to last excited state level of H 



binding energy of the second to last level is actually in- 
creased for laser intensities up to 0.25 kW/cnr^. However, 
at higher intensities the large positive lightshift of the 
lower levels dominates, pushing the second to last level 
to smaller binding energies and finally into the scattering 
continuum. 

In conclusion, a sufficiently strong 1/ long-range be- 
havior of the excited state potential is required to ensure 
a negative lightshift and compression of the wavefunc- 
tions with increasing intensity to allow for molecule for- 
mation via an optical FR. 



IV. STRONTIUM MOLECULE FORMATION 
VIA SHORT-PULSE PHOTOASSOCIATION 

While molecule formation via an optical Feshbach res- 
onance relies on field-free and field-dressed eigenstates 
and a time-dependence is only invoked to adiabatically 
connect these eigenstates, short pulse photoassociation 
is based on quantum dynamics on the much shorter 
timescale set by vibrational motion in the excited state 
potentials. A pump-dump photoassociation scheme is 
employed [13, [27| where the photoassociation pump pulse 
excites a small part of the density of colliding atoms from 
the electronic ground to an electronically excited state, 
cf. Fig. [TJd. The short pulse duration gives rise to a fi- 
nite spectral bandwidth. Population can be transferred 
into all excited state vibrational levels which are resonant 
within the bandwidth of the pulse, and a wave-packet is 
formed. The binding energy and spectral width of the 
wave-packet are determined by the pump pulse detun- 
ing and bandwidth [27j . These pulse parameters need to 
be chosen such that excitation of the atomic transition 
leading to strong loss of atoms from the trap ^33, 34] is 
minimized. This requires a large enough detuning and/or 
small enough spectral bandwidth. At the same time, the 
amount of population excited into molecular levels of the 
excited state, i.e. the photoassociation yield, needs to 
be maximized. The resonance condition together with 



the spectral bandwidth defines a 'photoassociation win- 
dow' around the Condon radius [IS, [1^, i.e. a range of 
internuclear distances where population is excited. The 
overall amount of density of colliding atoms within this 
windo w p oses an upper limit to the photoassociation 
yield (36[. At short internuclear distance corrsponding 
to large detunings, the population density becomes very 
small. A compromise leading to a large photoassociation 
yield while avoiding the excitation of the atomic tran- 
sition is found for detunings of a few wavenumbers. In 
this regime, pulse energies of a few nano-Joule were found 
sufiicient to completely photoassociate all available pop- 
ulation within the photoassociation window in the case 
of Rb2 and Cs2 molecules [271137}. 

Once the excited state wave-packet is created by the 
photoassociation pump pulse, it oscillates back and forth 
on the excited state potential. For heavy homonuclear 
molecules such as Rb2 or Sr2 and binding energies of a few 
wavenumbers, the vibrational periods are on the order of 
100 ps. Given excited state lifetimes of nanoseconds or 
more, the wave-packet undergoes many oscillations be- 
fore it decays due to spontaneous emission. A spatial av- 
erage of the wave-packet would hence be observed if the 
deexcitation to the electronic ground state proceeds by 
spontaneous decay, with partial wave-packets at short in- 
ternuclear distance yielding molecules and those at large 
internuclear distance producing pairs of hot atoms. In 
a time-dependent approach, a second short pulse is ap- 
plied when the wave-packet arrives at its inner turning 
point in order to transfer a maximum amount of popu- 
lation to the electronic ground state, cf. Fig. [T|3. The 
spectral bandwidth of the dump pulse can be chosen to 
populate a distribution of vibrational levels or a single 
ground state level [ij]. The efficiency of the dump step 
is crucially determined by the shape of the excited state 
potentials. Pure long-range potential wells with a softly 
repulsive part found in e.g. the 0~(P3/2) state of Cs2 or 
resonantly coupled excited states such as the 0+ states 
of Rb2 and Cs2 give rise to a pile-up of probability am- 
plitude at short internuclear distances and hence to an 
efficient dump step [13, . 



A. Dipole allowed transitions X^Eg 



The efficiency of photoassociation by the pump pulse 
relies on the 1 /R^ behavior of the excited state potential 
which is almost identical for Sr2 and Rb2. Similar pho- 
toassociation yields are therefore expected for Sr2 and 
Rb2. This is confirmed by Fig. [7] (see also Table I of 
Ref. [13): The final excited state populations of Sr2 are 
somewhat larger but of the same order of magnitude as 
those of Rb2 for similar pulse energies and detunings. 
The slight disparity in the final excited state populations 
is due to the difference in scattering lengths {as « 2 for 
®*Sr [m compared to as « 100 ao for *'^Rb) which re- 
sult in different ground state probability densities, i.e. 
in differing initial conditions, cf. Fig. [T]d. In order to 



FIG. 7: (Color online) Excited state population (filled cir- 
cles) and population in bound excited state molecular levels 
(open triangles) after the photoassociation pump pulse as a 
function of pulse energy for three different photoassociation 
pump pulse detunings ApA- The populations undergo Rabi 
oscillations as the pulse energy is increased. 
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time [ ps ] 

FIG. 8: Time-dependent projection of the excited state wave- 
packet onto ground state vibrational levels for a photoassoci- 
ation detuning of — 4cm~^ and a pump pulse energy of InJ. 
The upper panel shows the photoassociation pump pulse. 



compare the actual numbers of Pexc with those obtained 
for Rb2 [13 J an initial state with a scattering energy 
corresponding to 100 /xK was assumed in Fig. [T) Typi- 
cally, much lower temperatures for Sr2 are achieved by 
dual-stage cooling [2l| . Repeating the calculation for an 
initial state with a scattering energy corresponding to 
3^K, ApA = -4cm~^ and Eq = 2.5 x 10~^ a.u. yields 
a final excited state population which is reduced from 
1.8 X IQ-'^ in Fig. [3 to 9.8 x lO^^. Note that this does 
not necessarily imply a smaller number of molecules pro- 
duced by photoassociation, since the lower temperatures 
come with a higher density. The number of molecules per 
pulse depends on the specific configuration of the tra p, in 
particular on the density and the number of atoms [36l |. 
From the above considerations one can roughly estimate 
that at the least a similar number of molecules per pulse 
as for Rb2 (on the order of 10) can be expected for Sr2 

The pulse durations in Fig. [7] were chosen to be 
tfwhm = 10 ps for detunings of ApA = — 4cm~^ and 
ApA — — 5cm~^. The corresponding small spectral 
bandwidth of a lOps-pulse, Aw — 1.5 cm~^, ensures that 
no atomic transitions are excited. For /S.pA — — 19cm^^, 
the requirement of avoiding atomic transitions is also 
fulfilled by shorter pulses with larger spectral band- 
width. Since the spacing between vibrational levels in- 
creases with binding energy, a larger spectral bandwidth 
at larger detuning allows to address sufhciently many vi- 
brational levels to create a wave-packet. The pulse du- 
ration for ApA = — 19cm~^ was therefore taken to be 
TpwHM = 5 ps (corresponding to a spectral bandwidth of 
Aw 2.9cm-i). 

As the pulse energy is increased, Rabi oscillations in 
both the total excited state population, P'xc'*'' and the 
excited state population in bound levels, P^'^™'^'^, are ob- 
served in Fig. [T] At larger detuning, the free-bound 



Franck-Condon factors are significantly smaller. Since 
the Franck-Condon factors determine the Rabi frequency, 
ri(i) = fJ-E{t), the largest pulse energy of Fig. [7] corre- 
sponds to a 27r-pulse for ApA — — 19cm^^ compared to 
a 47r-pulse for ApA = — 4cm^^ and ApA = — 5cm~^. 
At small detuning (ApA = — 4cm^^), pulse energies of 
10 nJ and above lead to considerable power broadening. 
This is indicated in particular by the difference between 
the red circles and the red triangles which corresponds 
to the excitation of unbound atom pairs. It emphasizes 
once more the necessity to avoid any spectral overlap of 
the field with the atomic resonance. 

In a weak-field pump-dump approach, the maximum 
amount of population that can be transferred from the 
excited state wave-packet to bound ground state levels 
is achieved by a 7r-dump pulse [14| . It is determined 
completely by the (time-dependent) projections of the 
excited state wave-packet onto the ground state levels, 
PJ(t) = |((^^|^'e(t))|2. Figures [5] and m display the pro- 
jections onto the last least bound ground state levels for 
photoassociation detunings of — 4cm~^ and — 19cm~"'^. 
The oscillations of the projections mirror the oscillations 
of the wave-packet in the excited state potential. After 
one vibrational period for ApA = — 4cm~"'^ and after four 
vibrational periods for ApA = — 19cm~^ the oscillations 
show a beat structure which is due to the wave-packet dis- 
persion. The different vibrational periods in Figs. [5] and [5] 
are immediately rationalized in terms of the wave-packet 
binding energy, i.e. the photoassociation detuning. 

At small photoassociation detuning a non-zero projec- 
tion is observed for the last bound ground state level only 
(Fig. [5]). The wave-packet dynamics rationalizing this 
finding are those shown in Fig. [TJ): even when the wave- 
packet is at its inner turning point, it is peaked at a fairly 
large distance, near i? = 60ao, and it shows a large dis- 
persion. Clearly, such a wave-packet has not got a large 
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FIG. 9: Same as Fig. [8] but for a photoassociation detuning 
of — 19cm~^ and a a pump pulse energy of 3nJ. 



overlap with bound levels, and the maximum probability 
for transfer to the ground state amounts to about 0.6%. 
Choosing a larger photoassociation detuning, comes at 
the expense of a much smaller photoassociation excita- 
tion efRciency, cf. Fig. [T] The maximum projections 
seen in Fig. [5] correspond then to a transfer probability 
of about 15% to the last bound level and about 1% to the 
second to last level. Ground state levels with substantial 
binding energies cannot be accessed this way. 



B. Transitions near the intercombination line 

The efficiency of the photoassociation pump pulse 
for transitions near the intercombination line is limited 
by the quasi- long-range behavior of the excited 
state potentials. The corresponding free-bound Franck- 
Condon factors are several orders of magnitude smaller 
than for dipole allowed transitions. While this could in 
principle be compensated by employing amplified pulses 
with pulse energies of tens to hundreds of /iJ, a fur- 
ther obstacle is posed by the largely reduced density of 
excited state vibrational levels close to the dissociation 
limit. Even for strong pulses, the photoassociation yield 
is therefore expected to be significantly smaller than the 
yield predicted for dipole allowed transitions and weak 
pulses. 



V. CONCLUSIONS 

The perspective of forming Sr2 molecules in their elec- 
tronic ground state by adiabatic ramping over an opti- 
cal FR and by short-pulse photoassociation were investi- 
gated. In both schemes, molecules are formed coherently. 
This is unlike photoassociation with cw lasers followed by 
spontaneous emission 23]. Moreover, the molecule for- 
mation proceeds in a uni-directional manner for both the 



optical FR and short-pulse photoassociation. This dif- 
fers from two-photon photoassociation with cw lasers [2lj 
where the time-reversal symmetry is not broken, creat- 
ing a superposition of unbound atom pairs and molecules. 
While molecule formation via an optical FR is applica- 
ble in very tight traps such as deep optical lattices, short 
pulse photoassociation is best adapted to shallow trap- 
ping as found in a magneto-optical trap (MOT). 

Both schemes are promising when applied to opti- 
cally allowed transitions proceeding via the S^Sj ex- 
cited state, and similar molecule formation rates as in 
Rb2 [T^ . [27I [3^ are expected. This is rationalized in 
terms of the similar C3 coefficients for Sr2 and Rb2. The 
failure of both the optical FR and short-pulse photoas- 
sociation for transitions near the intercombination line 
is also attributed to the respective C3 coefficient. The 
1 jB? behavior of the relevant triplet state potentials is 
caused by spin-orbit coupling and was found to be rather 
weak [2^. As a result, optical FRs become open-channel 
dominated where the atoms are pushed apart rather than 
pulled toward each other, and short-pulse photoassocia- 
tion is inefficient due to lack of sufficiently many bound 
levels close to the dissociation limit. Molecule formation 
by photoassociation with cw lasers followed by sponta- 
neous emission might still be feasible. However, an esti- 
mate of the molecule formation rate based on accurate 
potential energy curves has yet to be given. 

While coherent molecule formation via dipole allowed 
transitions is predicted to be feasible, two notes of cau- 
tion must be made. (1) In the case of the optical FR, a 
very tight trap with corresponding short vibrational pe- 
riod is required to allow for an adiabatic ramp over the 
resonance while avoiding losses due to spontaneous emis- 
sion. This is due to the shorter excited state lifetime of 
Sr2 as compared to Rb2. (2) In pump-dump photoassoci- 
ation, the dump step is predicted to be significantly less 
efficient in Sr2 than in Rb2. The single, deep potential 
well of the Sr2 B^E^-state does not facilitate a pile-up of 
amplitude at short distance. This illustrates the difficulty 
of pump-dump photoassociation via a 'generic' potential 
where no mechanism for 'i?-transfer' such as resonant 
coupling between two electronically excited states is ef- 
fective. 

However, deeply bound levels could be reached in Sr2 
by a pump-dump scheme if additional coherent control 
techniques are employed. For example, resonant cou- 
pling can be engineered by applying an additional laser 
field. This has recently been suggested for short-pulse 
photoassociation in Ca2 |38|. The application of field- 
induced resonant coupling to strontium is beyond the 
scope of the present study. It is, however, expected to 
work for Ca2 and Sr2 alike due to their similar electronic 
structure. 

A second strategy to improve photoassociation and 
achieve a larger number of molecules aims at increas- 
ing the photoassociation yield, i.e. the efficiency of the 
pump step. While Feshbach-optimized photoassociation 
[39| is limited to the odd isotopes of strontium, flux 
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enhancement via coherent control of ultracold coUisions 
(40I l4ll . li^ could represent a general means of pushing 
the atom pairs closer together before photoassociating 
them. 

The molecule formation schemes presented here repre- 
sent an important first step toward deeply bound Sr2 
molecules which are required for future high-precision 
measurements Q . Once molecules in the last bound level 
are created, e.g. by ramping over an optical FR, they can 
be transferred toward deeply bound levels with a single 
strong, optimally shaped laser pulse HI, a train of weak, 
phase-locked pulses 1441 or stimulated Raman adiabatic 
passage (STIRAP) 
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